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r^ I Abstract 

\^ I In this paper, we study the dynamics and stability of a fundamental 

power system model when a time delay is imposed on the excitation of 
the generator. It is observed that sustained oscillations can arise in an 
otherwise stable power system through a delay induced Andronov-Hopf 
bifurcation. Numerical simulations are conducted to explore the dynamics 
of the time delayed system after the bifurcation which indicate period 
, doublings culminating in a strange attractor. 

S ! 1 Introduction 

f^ I Complex nonlinear phenomena in power system models have been extensively 

\0 ' studied in the past. Arnold diffusion, [1] and chaotic motions, [2] have been 

^D I studied in power system models described by the swing equations. The classic 

period doubling bifurcation route to chaos has been studied in a three node 

power system model, [3] using simplified models for the generator [4] , [5] , [6] , [7] 

^ ■ and later, with detailed generator models in [8]. Hard-limit induced chaotic 

^ I behavior in a single machine infinite bus formulation was studied in [9] by ex- 

• ^ ■ tensive numerical simulations where it was argued that the interaction of hard 

rN I limits and system transients led to sustained chaotic behavior. In [10], hard lim- 

jrt ' its were approximated by a smooth function to study bifurcations in the three 

node power system as an extension of [8] . 

While considerable attention ([1]-[10]) has been devoted to the study of power 
system models in the absence of time delays, not as much attention has been 
paid to understand power system dynamics in the presence of time delays with 
the exception of [11]. In modern power systems with digital control schemes, 
the acquisition, filtering and processing of signals can lead to possible time de- 
lays, [11]. In [11], the authors demonstrated an effective control technique using 
prediction to compensate for the delay in generator excitation. In the dereg- 
ulated power industry, communication delays arising from the use of phasor 



measurement units in wide area measurement systems could contribute to time 
delays in the acquisition of control signals. Thus, the focus of this paper is 
to study the effect of a time delay on the dynamics of a representative power 
system model. 

A Single Machine Infinite Bus (SMIB) system consisting of a generator inter- 
faced to the grid through a transmission line, Figure ^ is considered in this 
study. The SMIB system ( [12], [13]) has served as a long standing paradigm to 
understand the dynamics of a single generator when connected to a large power 
system. The simplicity of the SMIB system facilitates analysis and thereby 
enables one to gain insight in to the dynamic behavior of a synchronous gen- 
erator. The electro- mechanical equations of the generator are described by the 
standard flux-decay model and the excitation is represented by a single time con- 
stant system. Figure El The excitation voltage Efd is subject to a time delay 
T. The dynamics of the resultant system is then described by delay differential 
equations (DDEs). The reader is referred to [14], [15] for the general theory 
of DDE systems. In power system models described by ODEs, Hopf bifurca- 
tions usually arise when the loading of the generator, or control gains exceed 
certain critical values, [13]. However as remarked earlier, little attention (with 
the exception of [11]) has been paid to study power system models with delays. 
Therefore, we focus on understanding the role of a time delay on the dynamics 
of a synchronous generator. Later on, we show that even if the control gain and 
generator loadings are chosen to obtain a stable operating point, a critical delay 
in the excitation can cause the normally stable equilibrium to become unstable. 
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Figure 1: Single Machine Infinite Bus System 
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Figure 2: Excitation Control System 

2 System Model 

The system shown in Figure ^ illustrates a standard single machine infinite bus 
(SMIB) system equipped with a fast acting static excitation system (indicated 
by 'Exc'). Assuming a single time constant model for the excitation system 
shown in Figure |21 and the flux decay model for the generator, the dynamics of 
the single machine system in the presence of a time delay in excitation (Efd) 
can be described by, 
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(O, describes a limiter which imposes upper {E'Tf^) and lower {E^d") limits 
on the excitation voltage Efd of the synchronous machine. While the effect of 
excitation limits on the system dynamics and bifurcations has been previously 
studied in [9] and [10], the aim of the present paper is to understand the role 
of the delay term r. Therefore, in the rest of the analysis which demands 
smoothness of the models, the excitation limits are ignored which implies that 
Efdx — Efd- The rest of the variables {vd, id,Vq,iq) are given by, 
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The electrical power P^ and terminal voltage Vt are expressed in terms of the 
state variables from, 



Pe = Vdid + Vqiq, Vf = ^J vj + V, 



(8) 



3 Local Stability 

The local asymptotic stability of the equilibrium point for the system can be 
studied from the roots of the characteristic polynomial corresponding to the 
linearisation of S,- given by 
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In the linearisation of a single machine system expressed in (|1U|) , the parameters 
Ki, . . . Kq are known as Heffron-Phillips constants, [12], [13] which depend on 
the equilibrium and the machine parameters (see Appendix A). Substituting for 
A from I|1U|) in to yields the characteristic polynomial for the system which 
can be expressed as 



/,(A) = /i(A)+e-^V(A)=0 



(11) 



where. 



A (A) = A* + A3/3i + A2/32 + A/33 + (3^ 
gi{X)=X^p5+Pe 



The coefficients (3i,. . . f3e are given by. 
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with 

K[ = K,/2H, K'^ = K2/2H, K', = I/K^t'^,, 
K'^ = Ki/T'^^, K', = KaK^/Ta. K'e - KaKs/Ta 

In the absence of delay, by setting r = in (|ll|l we note that the characteristic 
polynomial is described by /o(A) = /i(A) + .gi(A) = 0. In the present study, 
the control gain, machine loading and other parameters are chosen so that the 
equilibrium point of the system without delay is locally asymptotically stable. 
With this assumption, the coefficients Pi,i — 1 ... 6 are all positive for typical 
values of machine parameters and loading conditions. In the following section, 
we shall study the local stability of the equilibrium as the delay r is varied. 

4 Delay Induced Oscillations 

Local stability of the equilibrium point of St- is governed by the roots of the 
characteristic polynomial /r(A) = 0, (|ll|l . Continuous dependence of the roots 
of Hll|l on T implies that there exists Tc such that 3?e{A} < when r G [0, Tc). A 
loss of asymptotic stability of the equilibrium occurs when Ke{A} = for some 
critical value of r = Tc > 0. Suppose the characteristic polynomial ifTT)) has a 
pair of purely imaginary roots A = ±iuJc when t = Tc^ we get 

/iM + e-'"^9iM = (20) 

Substituting for /i(.) and gi{.) from ((T^ and ((T^ we get, 

/3c + f3a cos(wt) = (21) 

Pb - Pa sin(cc;r) = (22) 

where, 
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From H21|l and H22|) . the pair {lOctTc) for which the system (|11|) has a pair of 
purely imaginary roots can be found from 
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By implicit differentiation of pi|l . we get 



dr /{(A)+gi(A)-re-^-5i(A) 
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Setting A = iojc and r = Tc in (|28|l . using H2()|) . (|27|l and simplifying the terms 
yields 



..S.{^)|_ = ?^fA§^ft) (29) 



where, 



D^ = Dl + Dl (30) 

-Da = /33 - 3/3iW^ - PaTc COs(WcTc) + 2Wc/32 sin(wcTc) (31) 

Db = 2tJc/32 - 4wj? + /3a Tc sin(cJcT) + 2ijJcl32 COs{ljJcTc) (32) 

Therefore, the sign of u is governed by the terms in the numerator of (j29(l . 
Substituting for the terms (32, Pb from Appendix B we get, 

^^-''^^ T-^K.T^ ^^nr (^^) 

Under normal operating conditions for a power system, the Hcffron-Phillips 
constants satisfy, 

Ki>0 (34) 

< i^a < 1 (35) 

< i^e < 1 (36) 

In addition, fast acting excitation systems typically have a small time constant 
so that < Te < 1. In view of the above conditions, we note that 

(/32 - /^s) > (37) 

The equations for /3a (wc) 7^ turn out to be complicated, making explicit analy- 
sis difficult. However numerical solutions of 1)26(1 over realistic ranges of operat- 
ing conditions yield nonzero values for the term Pa at the critical point in which 
case, cr > 0. Therefore, the results of the analysis so far can be summarized in 
the following proposition. 

Proposition 1 : The system S,- satisfies the sufficient conditions for a delay 
induced Andronov-Hopf-bifurcation under the following assumptions. 

[Al ] the system without delay Sq is stable ^ fo{^) — has no roots on the 
right half plane. 



[A2 ] if r 7^ 0, then A = is not a root of /^(A) l(TT|) . 

The assumptions [Al] and [A2] indicate that at the critical point, the character- 
istic polynomial has a purely imaginary pair of roots and all other roots of the 
polynomial have strictly negative real parts. The delay induced loss of operating 
point stability leads to the emergence of periodic solutions. The stability of the 
resulting periodic motion can be inferred from the sign of a certain quantity 
which proves to be too difficult in this case. Therefore, we resort to numerical 
simulations in the following section to study the dynamics of the system after 
delay induced instability. 

5 Simulation Results 

All numerical integrations were carried out using dde23 in MATLAB. A numer- 
ical simulation of the system E^ is shown in Figure 13 assuming typical values 
for the machine loading and control gain Ka (supplied in Appendix A) that 
results in a stable operating point in the absence of delay. Then the delay r is 
increased gradually from zero up to the critical value Tc when the equilibrium 
point loses stability. The rotor angle 5 is plotted in all time domain simulations. 
The transition from damped oscillations (when < t < Tc) to sustained oscil- 
lations (when T > Tc) can be noted from Figure |21 Incidentally, the solution of 
(|26|) and (|27|l yield Tc — 0.0533 which closely matches the critical delay observed 
by simulation, i.e. Tcri = 0.055. When the delay r is increased further, a period 
doubling (upper panel of Figure 0)) is noticed which eventually culminates in a 
strange attractor (lower panel of Figure 2)) . 
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Figure 3: The onset of oscillatory activity in the time delayed SMIB system. 
The critical delay Tc — 0.055. Upper panel shows a stable trajectory when 
T = 0.054 < Tc- Lower panel shows sustained oscillations when t = Tc- 
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Figure 4: Simulation of system trajectory beyond Tc- Upper panel indicates a 
period doubling when t = 0.059 followed by a strange attractor in the lower 
panel when r — 0.061. 



6 Conclusions 

Power system models are capable of displaying a rich variety of nonlinear phe- 
nomena which have been formally studied by bifurcation theory. In particular, 
the emergence of oscillations in power system models has been attributed to 
Hopf bifurcations which have been shown to arise when the control gains or 
machine loading exceed critical values. In this paper, the dynamics of a funda- 
mental power system model is studied when the field voltage of the generator 
is subject to a time delay. It is shown by preliminary analysis that a nominally 
stable operating point can be destabilized when the time delay exceeds a criti- 
cal value Tc through a Andronov-Hopf bifurcation. The behavior of the system 
for r > Tc is further studied by numerical simulations which indicate that the 
system undergoes a period doubling which eventually culminates in a strange 
attractor. Further analysis to determine the bifurcations of the ensuing limit 
cycle would be an interesting subject for future research. 
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Appendix A 

The Heffron-Phillips constants for a lossless network are given by 

^ ^ Ego COs((5o) ^ {Xq ~ x'^)iqo SVLl{5o) 
\Xq -t~ Xqj Xe r X^ 

Xe ~r X^ X(i ~\- Xq 

K, = (^^-^'<i)^';-iSo) (40) 

Xe n~ X^ 

_ -XqVdoCOs{So) _ -x'^Vqo sin((5o) 

{Xe + Xq)Vto {Xe + x'^)Vto 

where the subscript 'o' indicates the values of the variables at the equilibrium 
point which can be solved by setting (^- Q) to zero. 



Appendix B 

System Parameters : 

Xi = 1.79, xq = 1.71, x'^ - 0.169, P^ - 0.9, T^„ - 4.3, H = 5, wj, 

377, Xe = 0.4, i?e = 0. 

Exciter Parameters : 
Ka = 200, Ta = 0.05. 
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